Conditional Spontaneous Spin Polarization and Bifurcation 
Delay in Coupled Two- Component Bose-Einstein Condensates 



Chaohong Lee 1 , Wenhua Hai 2 , Lei Shi 1 , and Kelin Gao 1 
1 Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of 
Physics and Mathematics, The Chinese Academy of Sciences, Wuhan, 430071, People's Republic 

of China 

2 Department of Physics, Hunan Normal University, Changsha, 410081, People's Republic of 

China 
(February 1, 2008) 



Abstract 

The spontaneous spin polarization and bifurcation delay in two-component 
Bose-Einstein condensates coupled with Raman pulses are investigated. We 
find that the bifurcation and the spontaneous spin polarization depend not 
only on the system parameters, but also on the relative phase between two 
components. Through bifurcations, the system enters into the spontaneous 
spin polarization regime from the Rabi regime. We also find that bifurcation 
delay appears when the parameter is swept through the static bifurcation 
point. This bifurcation delay is responsible for metastability leading to hys- 
teresis. The area enclosed in the hysteresis loop increases with the sweeping 
rate of the parameter. 
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Electronic and nuclear -spin polarization in an atomic vapor with optical pumping have 
been investigated extensively [|TJ. Under conditions in which electronic spin exchange takes 
place faster than spin relaxation, spontaneous spin polarization appears. This interesting 
phenomenon is very similar to ferromagnetism and has been observed in wide ranges of 
atomic intensity, pump laser frequency and intensity. The appearance of spontaneous spin 
polarization means that the atomic vapor has two stable states with large spin polarization. 
The experimental realization of it has been applied to the field of optical bistability |§ . The 
atomic spin polarization exhibits striking hysteresis in switching between the bistable states 
0. This is analogs to those ferromagnetic system displays magnetic hysteresis 0. 

All the previous works have only considered the case of the thermal atoms, the exper- 
imental realization of multi-component Bose-Einstein condensates (BECs) 0J^] causes our 
interest in considering the similar behavior of the ultracold atoms. There are many differ- 
ences between thermal and cold atoms. The first one is that the collision among thermal 
atoms is noncoherent. However, when the temperature is close to the critical temperature 
realizing Bose-Einstein condensate (BEC) (T ~ Tbec)-, the collision among ultracold atoms 
is coherent due to the path between such collision is smaller than the phase coherence length 
|||7]]. Does this coherent property play a role in the polarization process of the ultracold 
atoms? Another difference between thermal atoms and cold ones is that the interaction 
strength of cold atoms can be controlled easily f|||. In this letter, we firstly show the co- 
herence among ultracold atoms gives rise to the conditional spontaneous spin polarization 
which is determined by both the system parameters and the relative phase. Then, we present 
the bifurcation induced by tuning the s-wave interaction between ultracold atoms. Lastly, 
we propose an experiment to confirm our prediction. 

We consider that the same kind of bosonic atoms, which are trapped in a single-well 
potential, are condensed in two different hyperfine levels |1 > and \2 >. Raman transitions 
between two hyperfine states are induced by the laser fields with the effective Rabi-frequency 
Q and a finite detuning 5. The internal Josephson effects PJ5|J^,pi]|, pT] , P^ | , coherent coupling 



effects [13 1, vortices [14| and spin waves [HI in such systems have stimulated great interest 



of many theoretists and experimentists. In the rotating frame, neglecting the damping and 
the finite-temperature effects, the coupled two-component BECs system can be described 



by a pair of coupled GPEs flO|,[T3|,|T 



ih d*i(r£ = ( H o + h mf + S)^, t) + f * 2 (r , t). 

Here, the free evolution Hamiltonians Hf = + V^(r) (i = 1,2) and the mean- field 

interaction Hamiltonians H^ IF = ^(a li \^ l (r , t)\ 2 + Oy|*j(r , t)\ 2 ) = 1,2, i ^ j) ( 
is the scattering length between states i and j which satisfies aij = ciji ). For the 
cigar-shaped trap (the trap frequencies satisfying u x = oj y » u z ), when the coupling is 
weak enough (i.e., the Rabi frequency satisfies fl/u z <C 1), the macroscopic wavefunctions 
can be written as the variational ansatz \I/j(r,t) = -?/>j(t)$.j(r ) (i = 1,2) with amplitudes 



V'i(^) = \J Ni(t)e toli ^ and spatial distributions $j(r). The symbols Ni(t) and a,(t) are the 
atomic population and phase of the i-th component, respectively. Due to the coupling is 
very small, the spatial distributions vary slowly in time and are very close to the adiabatic 
solutions to the time-independent uncoupled case for GP equations (1), being slaved by the 
populations [[K]]. Thus, the amplitudes obey the nonlinear two-mode dynamical equations 

HMt) = (e° 2 - 1 + u 22 \M)\ 2 + u 21 \M)\ 2 )M) + f ^(t), 



(2) 

ihf t ^(t) = {El + 1 + u u \M)\ 2 + u l2 \M)\ 2 )M) + f^M 

The parameters Ef = J $;(r )Hf^ i (r )d r, U {j = J \%(r )| 2 |$ j (r )\ 2 d r= U j{ and 

K — Q j <3>i(r )<3> 2 ( r )<^ r {hj — 1,2). The terms in K describe the internal tunnelling 
between two BEC states, whereas the terms in C/y, which depend on the numbers of atoms 
in each BEC state, describe the mean-field interaction between atoms. When U 2 \ and 5 
equals zero, these coupled equations can also describe the BECs in a double-well potential 
Introducing the Bloch's spin vectors 



_ + mi _ - Mi- Mi _ m? - (>,\ 

Obviously, u 2 + v 2 + w 2 = 1. When the total atomic numbers Nt = N\ + N 2 = + ip 2 *ip 2 
is conserved, setting the Planck constant h — 1, the Bloch's spin vectors satisfy 
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In this Bloch's equation, the parameters satisfy 7 = E\ — E\ + Nt(U22 — Uu)/2 — 5 and 
G = Nt{U%2 + U\\ — 2Ux2)/2. Comparing the above equation with the one for the linear 
case (Uij = 0) of equation (2), one can find that the mean-field interaction induces a shift 
(Gw) in the transition frequency and this shift is proportional to the relative population w. 

Taking |1 > as spin-up state and |2 > as spin-down state, the above two-component BECs 
system can be regarded as an ensemble of quantum spin- 1/2 particles. Thus, the longitu- 
dinal component w of the pseudospin describes the relative population, and the transverse 
components u and v characterize the coherence. In this language, the effective Rabi fre- 
quency causes an effective transverse magnetic field K along axis-w, the effective detuning 
induces an effective longitudinal magnetic field 7, and the mean-field interaction brings an 
effective longitudinal magnetic field Gw which depends on the longitudinal spin component. 

From the definition of the Bloch's spin vectors, we know that the above system can be 
described with only two independent variables. If we use the longitudinal spin component 
w and the relative phase = a<i — a.\ as independent variables, rescaling the time Kt to t, 
the motion equations 



dw/dt = — VT 



w 2 sin . 



d(J)/dt = -j/K - (G/K)w + w cos (j)/VT^ 



(5) 



are equivalent to the Bloch's equation. The above equations are consistent with those derived 
from the secondary quantized model ||. 

Below, let us analyze the stationary states of the coupled two-component BECs system. 
The stationary states can be obtained from the stable fixed points of the system. The fixed 
points correspond to those solutions satisfying dw/dt = and d(j)/dt = 0. In the region 
[0, 2tt) of the relative phase, we find two different modes of stationary states existing in the 
system: one is the equal-phase mode with zero relative phase (0 = 0), the other one is the 



anti-phase mode with n relative phase (0 = n). 

The number of the fixed points and the stationary states depend on the ratios 7/if, 
G/K and the relative phase. For the equal-phase mode, only a fixed point exists when 
G/K < 1 and this fixed point is stable. When G/K > 1, there are two stable fixed points 
and an unstable one for (G/K) 2 / 3 — (7/if) 2 / 3 > 1; there is only one stable fixed point 
for (G/K) 2 / 3 — (7/if) 2 / 3 < 1; the saddle-node bifurcations occur at the points satisfying 
(G/K) 2 / 3 — (7/if) 2 / 3 = 1. In the left column of Fig. 1, for the equal-phase mode, we show 
the values for the longitudinal component of the fixed points with different ratios 7/if and 
G/K. For the anti-phase mode, the corresponding case is different. When G/K > — 1, 
only a fixed point appears and it is stable. When G/K < —1, two stable fixed points and 
an unstable one exist for (G/K) 2 / 3 — (7/if) 2 / 3 > 1; only one stable fixed point emerges 
for (G/K) 2 / 3 — (7/if) 2 / 3 < 1; the saddle-node bifurcations occur at the points satisfying 
(G/K) 2 / 3 — (7/if) 2 / 3 = I. The fixed points of the anti-phase mode with different ratios 7/if 
and G/K are exhibited in the right column of Fig. I. In the Fig. 1, the fixed points between 
a pair of corresponding bifurcation points are unstable and the values for d(j/K)/dw at 
the bifurcation points equal zero. From the previous analysis, we find bistability exists in 
either the equal-phase mode or the anti-phase mode when the parameters obey (G/K) 2 / 3 — 
(7/if) 2 / 3 > 1. The appearance of bistability indicates the existence of spontaneous spin 
polarization (< w >= J T wdt/T ^ 0, T is the period for the oscillation of w) in this coupled 
two-component BECs system. When \K/G\ < 1, and 7/ if goes through the bifurcation 
points which satisfy (G/K) 2 / 3 - (7/if) 2/3 = 1, the spin polarization of either the equal- 
phase mode or the anti-phase mode is discontinuous at the bifurcation points. This means a 
first-order phase transition occurs. For the zero effective detuning 7, two metastable states 
with negative and positive spontaneous spin polarization coexist. As \K/G\ is increased 
to 1, the spontaneous spin polarization vanishes: this corresponds to a second-order phase 
transition. 

Similar to the thermal atoms, the spontaneous spin polarization can be induced by 
adjusting the coupling lasers. Additionally, the collisions among ultracold atoms can also be 
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controlled. Below, we will consider the bifurcation and spontaneous spin polarization in Bose 
condensed atoms induced by the ultracold collisions. Tuning the coupling laser with fixed 
intensity to a certain detuning satisfying 7 = 0, the bifurcation and the spontaneous spin 
polarization caused by the ultracold collisions in the above system can be obtained. For the 
equal-phase mode, only one stable fixed point w = exists when G/K < 1, two new stable 



fixed points w± = ±y 1 — {G/K)~ 2 appear and the original one w = becomes unstable 
when G/K > 1. This means a Hopf bifurcation occurs at G/K = 1. The system goes 
from the Rabi regime {G/K < 1) into the spontaneous spin polarization regime {G/K > 1) 
through the Hopf bifurcation. However, for the anti-phase mode, the Hopf bifurcation occurs 
at G/K = —1, there is only one stable fixed point w = for G/K > — 1, there two table 



fixed points w± = ±yl — {G/K)~ 2 and an unstable one w = for G/K < — 1. The Hopf 
bifurcations of both the equal-phase mode and the anti-phase mode are shown in Fig. 2. 
The solid lines are stable equilibria (stationary states), the dot lines are unstable equilibria. 

The Hopf bifurcations obtained from analyzing the equilibria in the previous are static 
bifurcations. Imagine now that the parameters are swept through the static bifurcation 
points. An interesting phenomenon emerges: the system starting close to the initially stable 
equilibrium does not immediately react to the bifurcation. Furthermore, it remains for some 
time close to the unstable equilibrium, then fast falls into one of the newly formed stable 
equilibria. This has been named as bifurcation delay which has been found in a variety of 
physical systems ||17|| . The bifurcation delay, which might lead to hysteresis, is the response 



to the bistability. For the equal-phase mode, fixed the effective detuning 7 = 0, slowly 
sweeping up the ratio G/K from R with sweeping rate r (i.e., G/K = R + rt, 1 ^> r > 0), 
choosing Rq < 1 and the initial state close to the equilibrium, the system evolves along 
the unstable equilibrium for a period of time after the ratio sweeping through the static 
bifurcation point {G/K = 1), then it quickly goes into a small oscillations around one of two 
new stable equilibria. The equilibrium, which the system evolves around lastly, determines 
by the state at the static bifurcation point. The system evolves around the up branch lastly 
when this state is close to the up branch; otherwise, the system evolves around the down 
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branch. When Rq > 1, slowly sweeping down the ratio through the static bifurcation point 
with initial state close to one of two stable equilibria, the system evolves near the stable 
equilibrium before it sweeps through the static bifurcation point, then it goes into a small 
Rabi oscillation around the ordinary equilibrium (w = 0). For the same sweeping rate, 
averaging the small oscillations, the process of sweeping up and down generates a hysteresis 
loop in the plane extended by G/K and w. The area enclosed in the hysteresis loop increases 
with the sweeping rate. This means that the energy exchanged between the atoms and the 
environments increases with the sweeping rate. The bifurcation delay in the equal-phase 
mode with different sweeping rate is shown in Fig. 3. For the anti-phase mode, a similar 
behavior can be observed near the static bifurcation point G/K = —1. 

To observe the spontaneous spin polarization and the bifurcation delay predicted in 
the above, we suggest a experiment based upon the works of JILA |5|]. Two BECs in the 
\F = 1, mj? = — 1 > and |2, 1 > spin states of 87 Rb are coupled by a two-photon pulse with the 
two-photon Rabi-frequency Q and a finite detuning 5. Thus, the control of the parameters 
K and 7 can be realized by adjusting the Rabi frequency and the detuning of the coupling 
lasers, respectively. The tuning of the parameter G can be accomplished with the Feshbach 
resonance ||. The longitudinal and transverse spin components can be measured with the 
state-selective absorption imagining and the Ramsey interference, respectively |15 . 



Summary, due to the coherent ultracold collision among condensed bosonic atoms, the 
bifurcation and the spontaneous spin polarization in coupled two-component BECs are de- 
termined by both the relative phase and the parameters. These phenomena are different 
from those only determined by the parameters, we name them as conditional bifurcation 
and conditional spontaneous spin polarization, respectively. For the zero effective detuning 
7, the Hopf bifurcation and bifurcation delay can be induced by the Feshbach resonance in 
either the equal-phase mode or the anti-phase mode. The system falls into the spontaneous 
spin polarization regime from the Rabi regime after the bifurcations occur. The appear- 
ance of bifurcation delay indicates the existence of metastability and hysteresis. Because 
of the inherently quantum coherence and superposition of two components, this quantum 



hysteresis might open the door to storage quantum data with Bose condensed atoms |L 

Acknowledgement 

The authors are grateful to the help of Dr. J. Feng in our institute. The 
work is supported by the National Natural Science Foundations of China (No. 
10275023, 19904013 and 19904014) and foundations of the Chinese Academy of 
Sciences. 

Figure caption 

Fig. 1 The fixed points for the system with different ratios "f/K and G/K. 
The numbers labelled on the lines are values for G/K. 

Fig. 2 The static Hopf bifurcation and the spontaneous spin polarization. 

Fig. 3 The bifurcation delay in the equal-phase mode for different values of 
sweeping rate which are labelled on the lines. 
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